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Abstract. Upper bounds are derived for codes in Stiefel and Grass- 
mann manifolds with given minimum chordal distance. They stem from 
upper bounds for codes in the product of unit spheres and projective 
spaces. The new bounds are asymptotically better than the previously 
known ones. 



1. Introduction 

Use of multiple transmit and receive antennas essentially increases the 
spectral efficiency of wireless systems (see and references therein). Analy- 
sis of Rayleigh flat-fading multiple-input multiple-output (MIMO) scenarios 
with m transmit antennas and n transmitted symbols, reveals that relevant 
coding schemes can be designed as collections of elements (points) in the 
complex Grassmann manifold - the set of m-dimensional linear subspaces 
in C n , if the channel is unknown to the receiver, and in the complex Stiefel 
manifold - the set of m orthonormal vectors in C n , if the channel is known 
to the receiver. An appropriately defined distance measure between the 
points characterizes diversity of the designed scheme. Following standard 
for coding theory considerations, we study the relation between the number 
of points (the size of a code) and the minimum distance between distinct 
code points. Our aim in this paper is to obtain new upper bounds for the 
size of codes in Grassmann and Stiefel manifolds. 

The most powerful technique for this kind of problems is the linear pro- 
gramming method (called also the polynomial method), initiated by Delsarte 
[2] . The method is very well understood in the case of 2-point homogeneous 
spaces (defined in the next section), where very explicit bounds, and also 
good asymptotic bounds on the rate of codes have been derived. Examples 
are the Hamming and Johnson schemes, treated in the unit sphere of 
W 1 [I], and the projective spaces [I], [Sj. 

When the underlying space is homogeneous and symmetric but not 2- 
point homogeneous, the situation is much more complicated, although the 
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principles of the linear programming method remain valid. The difficulties 
come from the fact that the zonal functions defined for these spaces are not 
functions of one variable, but afford several variables. The Grassmann spaces 
considered in this paper fall into this category. An attempt to overcome 
this problem was carried out in j^j. An asymptotic bound for the rate of 
Grassmannian codes was obtained, involving the asymptotics of the largest 
eigenvalue of some symmetric endomorphism. This bound however is not 
optimal since it was improved for m > 1 by some volume-type arguments 
for a large range of values of the minimal distance \7 a . 

There is one trivial case of symmetric spaces of rank m > 1 for which the 
classical treatment of the linear programming method is easily extended: 
it is the direct product of 2-point homogeneous spaces, such as the direct 
product of m copies of the unit sphere. An example of a similarly easy case 
is provided by the non-binary Johnson space [SJ 03 , that is the product 
of the Hamming and the binary Johnson scheme. 

The approach developed this paper is to relate Grassmann and Stiefel 
spaces and their associated codes to various products of 2-point homoge- 
neous spaces, and hence to derive upper bounds for these codes in a rather 
easy way. The asymptotic versions of the new bounds provide the best 
currently known asymptotic bounds. 

The paper is organized as follows. Definitions and known results are 
given in Section [21 Section |3] describes various relations between the spaces 
and their codes. The simplest one connects Grassmann and Stiefel spaces 
to the unit sphere of an asymptotically equal dimension; therefrom, for 
example, we obtain a bound for the asymptotic rate of Grassmannian codes 
that already improves upon the previous ones (Theorem IM.2J) . Section 0] 
develops the Delsarte polynomial method for the products of spaces under 
consideration, including the classical method that involves the Christoffel- 
Darboux formula, and derives upper bounds for the size of the associated 
codes. A bound for the asymptotic rate of these codes is obtained. Section 
El discusses the consequences for the Grassmannian and Stiefel codes. In 
particular, we show that the bound obtained on the asymptotic rate of 
Grassmannian codes from the product of projective spaces is sometimes 
better than the one obtained in Sectional We conclude in Section El 

2. Definitions and known results 

We shall use the following notations and definitions. We say that f(n) < 
g(n), f{n) ~ g(n), f(n) > g(n) if lim^oo ^ < 1, lim^oo ^ = 

1, linin^oo > 1, respectively. A code in a metric space (X, d) is a finite 
set contained in the space, and a codeword is an element of the code. The 
size of a code C is its cardinality, and its rate is R(C) := ^ In \ C\ where In 
denotes the natural logarithm. The meaning of n will be defined separately 
for each space. Indeed, a more consistent and general definition of the rate 
of a code in a manifold X would be di ^ X ) m 1^1 ■ ^he minimum distance of 
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a code is the minimum distance (induced by the relevant metric) between a 
pair of distinct codewords. A metric space (X, d) is called 2-point homoge- 
neous, if X affords the transitive action of a group G, such that the orbits 
of the action of G on X x X are characterized by the distance d. In other 
words, for all (x,y) E X and g E G, d(g(x), g{y)) = d(x,y), and moreover, 
for all pairs (x,y), (x',y') E G, there exists g E G such that g(x) = x' and 
g{y) = y' if and only if d(x,y) = d(x',y'). It is a well-known fact that the 
compact Riemannian manifolds that are two-point homogeneous are exactly: 
the unit sphere S" 1-1 , the projective spaces P Ti_1 (i ; C) where K = R, C, EI and 
the projective plane over the octonions P 2 (0) (see jllj . and ^2] f° r more 
about the octonions and P 2 ( 



2.1. The real compact two-point homogeneous spaces. The unit 
sphere of the Euclidean space P n is denoted S n ~ l , namely, 

(1) 5™- 1 :=|(x 1 ,...,x n )ElR"| X>? = lj- 

The standard scalar product in P n , given by (u-v) = Xa=i ^jUj, defines the 
Euclidean distance between two points of S n : 



(2) ll» 



II 



1=1 



The angular distance between u and v is defined by the angle 0(u, v) E [0, tt], 
also denoted when the context is clear. We have of course 

(3) \W-v\\ = v^Vl-costf. 

The best known asymptotic bounds on the rate of spherical codes as 
a function of the minimum distance are given in the following theorems. 
Theorem 12.11 is the Chabauty-Shannon-Wyner, which is the analog of the 
Gilbert- Varshamov bound for codes over finite fields, Theorem 12.21 is the 
linear programming bound, and Theorem 12.31 is an Elias-Bassalygo type 
improvement due to Yaglom. See ^3 El for details on these bounds. 

Theorem 2.1 (|I5j.[IB].p7j). Th ere exists a sequence of codes {C n }^L-L 
such that C n is of length n, minimum angular distance 9 < tt/2, and rate 
R(C n ) := -ln|C n | which satisfies 

(4) R(C n )>- In sin0. 

Theorem 2.2 (H]). Let C be a spherical code with minimum angular dis- 
tance 9 < tt/2. Then, when n — > oo, 

(5) «(C)<W):=ii# H A -™ 



2 sin 9 \ 1 + sin I 

where H is the entropy function, H{x) := — xlnx — (1 — x) ln(l — x) 
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Theorem 2.3 ([I]). Let C be a spherical code with minimum angular dis- 
tance 9 < tt/2. Then, for any 9/2 < <p < tt/2, when n — > oo, 

(6) R{C) < R(a) -In sin <p, 

where sin(a/2) = S1 °^^ ; and R(a) is an upper bound on the asymptotic 
rate of spherical codes with minimum angular distance a. 

For R = Rlp, the optimal choice of ip is discussed in [1]; it corresponds 
to a ~ 63° and gives a better bound than Rlp{9) when 9 is smaller than a: 



(7) R(C) < R Y (0) ■= - In Vl - cos 9 - 0.0686. 

We denote by Rs(9) the function that provides the best known bound for 
the asymptotic rate of spherical codes: 

' R Y {9) if 9 < a 
Rlp{9) ifa<9<ir/2 



(8) R S (9) 



Remark 2.4. The asymptotic rate of spherical codes with minimum angular 
distance at least tt/2 is known to be equal to zero. This is a consequence of 
the Rankin bound C|18|. see also |13| or |14j ). 

The other real compact manifolds which are two-point homogeneous can 
be treated in a similar way. As was recalled before, these are the projective 
spaces P n_1 (i ; C) where K = R, C,H (the field of real quaternions) and n > 3, 
and the projective plane over the octonions P 2 (0). In order to treat the 
fields of coefficients in a uniform way, we extend the definition of (x • y) so 
that, for all x,y E K n , (x ■ y) = Ya=1 x iVi^ where the conjugation x — > x 
is the standard one over K = C, H, O and is the identity over M. Also we 
conventionally assume that n = 3 when K = O. The group G under which 
these spaces are two-point homogeneous is respectively the orthogonal group 
0(M n ), the unitary groups U(K n ) with K = C,H, and the Lie group (see 
|12j for this last case). 

The angular distance between p and q in P n ~ 1 (K) is defined by the angle 
9 = 9{p, q) £ [0, tt/2] such that cos 9(p, q) = |(e • /)| where e, / are arbitrary 
chosen unit vectors of the lines p, q. It is shown in [3] and that the 
linear programming method applies to these spaces. The derived asymptotic 
bound for the rate can also be obtained from the bounds for spherical codes, 
because to a code C in ¥ n ~ l (K) one can associate a code in S cn_1 with the 
same size and a minimum angular distance at least equal to the one of C, 
selecting a unit vector in each element of C. One obtains: 

Theorem 2.5 ([!]). Let C be a code in F 1-1 ^) with minimum angular 
distance 9 < tt/2. Let c := 1,2,4 respectively when K = R, C,H (so that 
c = [K : R]J. Then, when n — > oo, 

(9) R(C):=-\n\C\<cR s (9) 

n 
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2.2. The Grassmann space. Let K be the real or the complex field. The 
Grassmann space Q m ,n{K) is the set of all subspaces of dimension m in 
K n . It is a homogeneous space under the action of either the orthogonal 
group 0(M n ) or the unitary group U(C n ). We will denote Q m ,n when K is 
arbitrary. It is worth noticing that when m = 1 we recover the projective 
space. Several metrics have been defined in Q m ,n, see ^3120]. In this paper 
we consider the chordal distance, which was introduced in [2^ and studied 
in |19| 1211 120| . The following two definitions for the chordal distance 
d c (p, q) are equivalent. 

Definition 2.6. Given the planes p,q E G m ,n, o,pply the following procedure. 
Initialize the sets of unit vectors A = and B = 0. In the ith step, choose 
the vectors cij , bi such that: 

(i) ai is contained in p and bi is contained in q. 

(ii) ai is orthogonal to all the vectors in A and bi is orthogonal to all 
the vectors in B. 

(iii) Among all the vectors satisfying the conditions in (i) and (ii), the 
angle between ai and bi is minimal (i.e., their inner / Hermitian 
product module is maximal). 

Set Oi to be the angle between ai and bi, insert to A and bi to B, and 
proceed to the next step, until m angles < 0± < 62 < . . . < 6 m < tt/2, 
called the principal angles between p and q, have been defined. Then the 
chordal distance is 



ni 

m — cos2 @i- 

i=l 

Lemma 2.7 (|21j). For a plane p £ Gm,n> ^ A p be a p x n matrix whose 
rows form an orthonormal basis of p, and let tt p := A*A P be the matrix of 
the orthogonal projection on p (A* denotes the Hermitian conjugate of A p ). 
Then, the projection matrix ir p does not depend on the choice of A p , and, 
given two planes p,q £ Q m ,n, the chordal distance is 



(10) d c (p, q) ■= y m — trace (ir p o Tr q ). 

We review some known bounds on the size of codes in Grassmann spaces. 
We recall that c = 1 if K = M and c = 2 if K = C. The first work in 
the area is due to Conway et al. They found an isometric embedding 

from G m ,nW to the Euclidean sphere with radius \J m{n ~ m) in mK™- 1 )(«+ 2 ) 
(a similar embedding exists also for {5 m ,n(<C), see ^2j)- This enabled them 
to use the Rankin bound on spherical codes |18j to derive the following 
non-asymptotic bound. 



d c (p, q) 



\ 



m 

E 

i=l 



sin 



\ 
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Theorem 2.8. Let C be a code in C/ mjn (R) with minimum chordal distance 
d. Then 

{m(n—m) |C| 
,/|C|>In(n + l) 

Later, an asymptotic expression to the volume of a ball in Grassmann 
spaces was derived by Barg and Nogin ^5], yielding the analogue of Gilbert- 
Varshamov and the Hamming asymptotic bounds for the rate R(C) := 
— In |C| of Grassmannian codes. 

Theorem 2.9 (19 ). For any pair of constants d,m, such that d < \fm, 
there exists an infinite sequence of codes {C n } in Qra n with minimum chordal 
distance d and rate 

(12) R > -cm In - 6 



"m 

Theorem 2.10 (QU). Let C be in Qrn.n with minimum chordal distance d. 
Then, when n — > oo, 

(13) R(C) < -cm In (yj 1 - y/T Z d?/2m^j . 

A linear programming bound was derived by Bachoc [H]. 

Theorem 2.11 (0). Let C be a code in Qm,n{^) with minimum chordal 
distance d. Then, when n — > oo, 

(14) R{C) < m[(l + p) ln(l + p) - p In p], 

where 

(15) p = —m(y/m/d — 1). 

We note that the derivation of the linear programming bound is not 
a straightforward analogy to the derivation of this bound in other metric 
spaces, since it involves a family of orthogonal generalized Jacobi polyno- 
mials with several variables, and that the bound Q14|) in the case m = 1 
coincides with the bound for the real projective space of Theorem 12.51 The 
case of the complex Grassmann space is not treated in 0| but could be 
treated in a similar way. 

Finally, a new upper bound was introduced recently by Barg and Nogin 
[7j, using Blichfeldt's density method. 

Theorem 2.12. Let C be a code over the Grassmann space Gmnffi) with 
minimum chordal distance d. Then 



(16) R(C)<-m In 



It is immediate to see that this bound improves upon (fH3)) and upon ((11)) 
for m > 1 and for a large range of values of d. Hence, until this paper, 
Theorems 12.111 and 12.121 provide the best asymptotic known upper bounds 
on the rate of codes in Q m ,n{ 
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2.3. The Stiefel manifold. The Stiefel manifold V mn (K) is the set of 
m-tuples of orthonormal vectors in K n , or equivalently 

V m , n (K) = {X£ M mxn (K) | XX* = Id m }, 

where Id m is the m x m identity matrix. The orthogonal group 0(M n ) 
if K = R, respectively the unitary group U(C n ) if K = C acts transi- 
tively on V mtn (K), and this space can be identified with the set of classes 
0(R n )/0(R n '- m ), respectively U(C n )/U(C n - m ). 
The distance considered in coding theory is 



d(X, Y) := \\X - Y || = v /trace((X-y)(X*-y*)). 

In other words, d(X, Y) is the Euclidean distance between X and Y, when 
X and Y are regarded as one-dimensional vectors of length mn. We refer 
the reader to [22] for a treatment of codes in Stiefel manifolds. 

3. More spaces and their interconnections 

The simplest of these connections relate Grassmann and Stiefel spaces 
to a single unit sphere, and allow to apply directly the known bounds for 
spherical codes to the Grassmannian and Stiefel codes. We start with them, 
then we introduce the products of spheres and projective spaces and their 
relations with Grassmann and Stiefel spaces. 

3.1. G m ,n and S cmn ~ 1 . We follow the notations and definitions of Section 
12.21 For all p,q £ Q mj7l (K), we set 

m 

a (p, q) '■= c °s 2 &i = trace(-7Tp o 7r ? ). 
We define a mapping 

in the following way. We select for all p £ Gm,n(K), an orthonormal basis 
(ei, . . . , e m ) of p whose elements belong to K n . With the usual identification 
of C and M x R through the mapping z = x + iy — > (x,y), we consider these 
elements in R cn . Then (5{p) is chosen to be the element of R cmn obtained by 
the concatenation of e±, . . . , e m , divided by ^frn. Obviously, (5{p) £ gcmn-i_ 
The new bounds for Grassmann spaces rely on the following lemma. 

Lemma 3.1. For allp,q £ Q mn (K), 



cos 9(P(p),f3( q ))<J^^-. 

V m 

Proof. Let (3(p) = e, obtained from an orthonormal basis (ei, e2, • • • , e m ) of 
p and [3(q) = e', obtained from an orthonormal basis (e[, . . . , e' m ) of q. We 
compute cr(p,q) = trace(-7r p o ir q ). Let A p , A q denote the m x n matrices 
whose rows are the basis elements ej, e- respectively. Then 



iPi q) = trace(7T p o ir q ) = trace(A* A p A* A q ) = tiace(A p A*A q A. 
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The entries of the matrix A p A* are the hermitian products (ej • e' ). So 
we obtain: 

(17) a(p,q)= £ \(fii'4)\ 2 

l<i ,j<m 

If K = R, we obtain from Cauchy-Schwartz inequality 

(18) cos 9(e, e>) = (e • e') = < / SE53g < S 

m V m V m 

If = C, let us denote by 5ft(z) the real part of a complex number z. 
In the identification C n = M? n recalled above, the standard scalar product 
on M. 2n is given by y)). With the obvious inequality 5Ji(/i(x, y)) 2 < 

\h{x,y)\ 2 , we obtain the same inequality cos#(e,e') < \J a (where e and 
e' are considered in the unit sphere of R 2n .) 

□ 

Let us recall Definition 12.61 of the chordal distance d c (p,q) in Grass- 
mann spaces. The definition involves the construction of orthonormal basis 
(»!,..., a m ) and (&i, . . . , b m ) for p and g respectively, such that the princi- 
ple angles 6>j satisfy cos6>j = |(a, • and cr{p,q) = ^^cos 2 ^. Lemma 
13.11 shows that if one chooses arbitrary orthonormal basis (ei,...,e m ) of 
p, [e'i, . . . , e' m ) of q, and defines an alternative set of "principal angles" 
d[,...,9' m by 6[ = arccos|(ej • e^)|, then a(p,q) > Ys1Li cos2 Thus, an 
upper bound on the chordal distance between p and q is obtained. 

It follows from Lemma 13.11 that the bounds for spherical codes (of the 
sphere S cmn ~ l ) can be applied to codes in Grassmann spaces. We obtain 
for the asymptotic rate: 

Theorem 3.2. Let C be a code in Qyn^Yi^K^) with TniniTnal chordal distance 
d = \fm — s } and let 9 = arccos(y / 5/m). Then, when n — » +oo ; 

(19) R(C) < cmR s (e) 

where Rs is defined in ©. 

Figure n compares (|19|) with some of the existing bounds, given in Section 
12.21 The bound (|16j) can also be obtained using the mapping 0, joined with 
the Rankin-Blichfeldt bound [7] . It is well known that the Rankin-Blichfeldt 
bound |18j is improved by the bound Rs(0) for all values of 8. 

Remark 3.3. In |21| . the authors introduce an isometric embedding of 
G m ,nO&) into a unit sphere, but the dimension of this sphere is much larger 
than the one of the Grassmann space ((n—l)(n+2)/2 instead of m{n — m)). 
Instead, the dimension of S cmn ~ 1 is asymptotically equivalent to the one of 
Qm,n(K), but our embedding is not isometric. 
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Figure 1. Upper bound on the asymptotic rate of real 
Grassmannian codes with minimum chordal distance d, m = 
3. From top to bottom: a linear programming bound (|14|) . a 
Blichfeldt-type bound Q16|>. and the new bound (fTT3|) 



3.2. V mn and 5 cmn - 1 . 



Lemma 3.4. Let X, Y € Vm,n(-?0, i$T = R, C. Let (ei,...e m ) denote the 
rows of X, respectively (e[, . . . e' m ) for the rows ofY. Then 



d(X,Y) = V2 



m 



Proof. We calculate 

||X - y|| 2 = trace((X - F)(X* - Y*)) 

= trace(XX* - XY* - YX* + FF* 
= 2m - 2K(trace(XY*)) 



since XX* = YY* = Id m and trace(Xy*) = trace(XY*) = trace(yX*). 
We conclude with 



trace(Xy*) = ^Tfa ■ e-). 



i=i 



□ 



Again with the identification of C n with M? n , we view V niTI .(C) as a sub- 
manifold of V m> 2n(^) endowed with the distance 



d{X,Y) = \\X-Y\\ = V2 



m 



t=i 
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We see that the obvious mapping: 

7 : V m>n (K) - S^- 1 



X i y 7(X) = -^=(e 1} . . . , e m ) 



is this time, up to a suitable scaling of the distances, an isometry. Hence the 
bounds for spherical codes also apply to V mjn (if), probably in a quite effi- 
cient way. Still, one constraint is not encoded in it: the fact that the vectors 
ei are pairwise orthogonal and of norm 1. We resume these observations in 
the following theorem. 

Theorem 3.5. Any upper bound on the size of codes in S cmn ~ l with mini- 
mum angular distance 9 = arccos(s/m) is also an upper bound for codes in 
the Stiefel space V m ,n(K) with minimum distance \[2 V m — s. In particular, 
we have for the asymptotic rate, when n — > +oo, 

(20) R(C) < cmR s (8) 

where Rs is defined in (|SJ). 

3.3. Qm,ni Vm,n and products of spaces. So far we have established a 
relation between codes in Qran and Vrnn and codes in Q crnn i # It is worth 
noticing that the mappings j3 and 7 defined above factor out by (5 cn_1 ) m , 
since the elements are unit vectors. Hence bounds for codes in (S' cn ~ 1 ) m , 
will imply bounds for codes in Qra n and V772 This is the motivation to 
the generalization of the linear programming method to the product of unit 
spheres, and more generally to the product of 2-point homogeneous spaces, 
which is proposed in the next section. As we shall see, the asymptotic bound 
for the rate of codes in (S cn ~ 1 \ is not better than for 5 cmn_1 ) hence doesn't 
improve on (|19j) and 1)20(1 . A better result is obtained for Grassmann spaces 
with the product of projective spaces. 

We now define more precisely the products of spaces and their associated 
distances that will be studied in the next section. We start with the product 
of m copies of the unit sphere of M": 

(5"- 1 ) m = {e = (e 1 ,...,e m )|e,G^- 1 }. 

We consider on the distance given by 



d(e,e') 



'2 m 



f 119 



Em 
, _ i=l COS 



1 
m 



where cos^j = (e, • e'j). We attach to a pair e, e' G (S n 1 ) m an angle 
9 = 9(e, e') G [0, vr] such that 

E™ 1 cos 0i 



(21) cos) 

m 
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and we call 9 the angular distance between e and e'. The angle 9 is also the 
angle between the vectors e/^Jm and e'/^/m, viewed as elements of S mn . 

We define, for the remaining 2-point homogeneous spaces recalled above, 
and without specifying the field K, 

^ r -Y = {p={p 1 ,...,p m )\p l e¥ n - 1 }. 

We attach to a pair p,p' £ (P n 1 ) m an angle 9 = 9(p,p') £ [0,7r/2] such 
that 

2 a ^T=\ cos2 & 
cos 9 = — — 

m 

where 9{ = 9(pi, qi) and we call 9 the angular distance between p and p' . We 
consider on (P™ 1 ) m the "chorda!" distance given by 



d(p,p') 



m 

ni - cos2 Qi - v ? " SUI ' 
i=l 



In order to derive bounds for codes in Grassmann spaces G m ,n(K), we 
shall make use of the mapping: 

v:9 m , n (K)^(¥ n -\K)) m 

defined in the following way: for all p £ Gm,n(K), we choose a m-tuple 
(pi, . . . ,Pm) of pairwise orthogonal lines of p. We set v(jp) = (pi, . . . ,p r , 
Because of the equation (j!7[) . we have similarly: 



J mj- 



cos 9{v{p), f(q)) < 

m 

hence the bounds for codes in (P n_1 (X)) m apply to codes in Grassmann 
spaces. 



4. Bounds for codes in the product of 2-point homogeneous 

SPACES 

In this section, X denotes one of the spaces S n ~ l , P n_1 (X) where K = 
M, C, H, or the projective plane over the octonions P 2 (0). We derive bounds 
for codes in X m with a given minimum distance, following Delsarte's linear 
programming method as performed in 0]. As a reference on orthogonal 
polynomials, we refer to |23j . 

4.1. Review of the necessary material on the harmonic analysis of 
the spaces X. We recall that, to each of these spaces is associated a family 
of orthogonal polynomials of one variable, which are the zonal polynomials 
relative to the action of the group G (see jl], 0, |21])- For X = 5 n_1 , these 
polynomials are the Gegenbauer polynomials with parameter re/2 — 1 and 
associated orthogonal measure (1 — x 2 )(™~ 3 )/ 2 on the interval [—1,1]. For 
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X = P n 1 {K), these polynomials are Jacobi polynomials with parameters 
(a, (3) defined by: 

a = |(n-l)-l, = l~l. 
More precisely, the values of (a, (3) are as follows: 





c 


a 





R 


1 


(n-3)/2 


-1/2 


C 


2 


n - 2 





m 


4 


2n-3 


1 


o 


8 


7 


3 



The orthogonal measure associated to the parameters (a, (3) is x@(l — x) a 
over the interval [0, 1]. We generically denote by Pk{x) these polynomials, 
with deg(-Pfc) = k and -Pfc(l) = 1. We let fj,(x) denote their normalized 
associated orthogonal measure and [P, Q] the corresponding scalar prod- 
uct on R[x] (so that [P,Q] = J P(x)Q(x)fx(x)dx and [1,1] = 1). More- 
over, we have [Pfc,P&] = d^ 1 where dk denotes the dimension of the ir- 
reducible representation of G associated to Pk (e.g. when X = S 11 ^ 1 , 
d k = dimHarm fc = ("+J- 1 ) - (^If))- 

The three-terms relation expresses xPk(x) as a linear combination of the 
polynomials Pf. 

xP k (x) = a k Pk+i{x) + b k Pk(x) + c k Pk-\{x) 

for some sequences of rational numbers (a^), (bk), (c^). It is enough for our 
purpose to know that (a&) is bounded when n and k tend to +cx> with n/k 
tending to a finite limit. For example, when X = S n ~ 1 , 



dk 



n-2 + k 
n-2 + 2k' 



For all (u,v) € X, we define 
t(u, v) = 



(u-v) = cos 6{u,v) ifX = S n - 1 
cos 2 6(u,v) if X = ¥ n ~ 1 (K). 

The zonal function on X associated to Pk is given by: 

(u,v) -> P k (t(u,v)). 

The so-called 'positivity property' related to these polynomials, and un- 
derlying the linear programming method in X, is the following: for all code 
C CX, and for all k > 0, 



ueCvec 



BOUNDS FOR CODES IN PRODUCTS OF SPACES 



13 



4.2. The linear programming method on X m . Now we consider the 
product spaces X m . The positivity property generalizes to the following: 

Lemma 4.1. Let C C X m . Let us denote elements of C by u = (u\, . . . ,u m ) 
with m G X. For all (ki,..., k m ) £ N m , 

m 

uec vec i=i 

Proof. This is the positivity property in the product space X m . The group 
G m acts transitively on X m ; the G m -irreducible components of L 2 (X m ) are 
the tensor products of the G-irreducible components of each L 2 {X) and the 
associated zonal functions are given by the polynomials in the m variables 

•£"1 j • • • j x m 

m 

l\PkM), (h,...,k m )en m 
i=i 

in the way: 

m 

(u,v) I ^ Y[P ki (t(Ui,Vi)). 

1=1 

□ 

Remark 4.2. In a sense, the polynomials TYiLi Pki(xi) are fake multivari- 
ate polynomials since the m variables are separated. As we shall see, for 
this reason it is much easier to deal with them, compared with the zonal 
polynomials for the Grassmann space (see [Hj^. 

The polynomials ni^=i Pki(, x i) generate the polynomial algebra 
C[x 1} ... , x m ] , and are orthogonal for the product measure 

in 

8=1 

with support [— l,l] m when X = S" 1-1 , respectively [0, l] m otherwise, and 
where A is chosen so that the total measure is equal to 1. The associated 
scalar product on R[xi, . . . , x m ] is denoted by [, ]. We take the following 
notations: a multi-index in N m is denoted by k = (ki, . . . , k m ) and we 
define for x = (xi,...,x m ), Pu{x) = Pk(x 1 , . . . , x m ) := n^=i p ki(.Xi), and 
dk := YliLi d-ki- Obviously we have, for all k and /, 

(22) [Pk,P L ]=hid k 1 . 
Moreover, we define 

m 

a(x) := y^Xj. 
i=i 

For any angle 9 we denote 

f t = cos^ ifX = S n - 1 
| t = cos 2 9 if X = P n - 1 (K) 
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Now we can formulate the usual associated linear programming bound: 
Proposition 4.3. Assume F £ M.[x\, . . . ,x m ] satisfies the conditions: 

(i) F = fkPk with fk > for all k, and fo>0 
k 

■ IfX = S n -\ F( Xl ,...,x m )<0 for all (x\, . . . , x m ) G [—1, l] m 
.... I such that cr(x) < mcosQ = rat 

I IfX^S n -\ F(n i m )<0 for all (xx, . . . , x m ) 6 [0, l] m 

such that <j(x) < mcos 2 9 = mt 

Then, any code C in X m with minimum angular distance 9 satisfies 

|C|<%^> 

JO 

Proof. We reproduce the standard argument. Let 

S := ^2 5^ ^(*C u i> ?; i)> • • ■ ,t(u mi v m )). 

u£Cv&C 

The pairs (u, v) with u = v contribute in this sum for |C|.F(1, . . . , 1). From 
condition (ii) and the assumption that for u ^ v G C, cos9(u,v) < cos 9, 
the other terms are non positive. Hence, S < |C|.F(1, . . . , 1). 
On the other hand, we have 

S = ^2fki ^2 p k{~t{ui,vi),. . . ,t(u m ,v m )). 

k u,v£C 

The term corresponding to k = = (0, . . . , 0) gives /o|C| 2 while the other 
terms are non-negative from the positivity property of the polynomials Pk 
(Lemma l4.1jl . Hence S > /o|C| 2 . The two inequalities lead to the announced 
bound. 

□ 



4.3. Examples of small degree. Let us work out the case of polynomials 
of small degree. 

(i) X = S"- 1 

(a) Degree 1: we take F = (x\ + - ■ ■+x m ) — mt. Since P\{x) = x, F 
satisfies the hypothesis of Proposition 14.31 if and only if t < 0. 
We obtain: 

(23) Ifcos0 = t<O, |C|<l-j. 
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(b) Degree 2: we take F = ((xi + • • • +x rn ) + m) ((xi + • • • + x m ) 
mt). We have 

F = {x\ H h x m ) 2 + m(l - t) (xi H h x m ) - m 2 t 



+ 2 XjXj + m(l — t) x i ~ m2 t 

i<j 

y^(xf - — ) + 2 XjX 7 - + m(l - t) Xj + — - m 2 
^— ' n n 



Since Pa(^) = (x 2 — V n )/(1 ~~ l/ n )> ^ satisfies the hypothesis 
of Proposition 14,31 if and only if ^ — m 2 t > 0. We obtain: 

(24) Ifcos0 = t<J-, |C7|< 2mn(1 - t} . 

mn 1 — rant 

The two bounds take the value 1 + mn at their crossing point, 
corresponding to t = —1/mn. . 
(ii) X = P n_1 (i^), Degree 1: we have, up to a multiplicative factor, 
Pi(x) = x - a +~p\ 2 = x- \. We take F = (x\-\ h x m ) - mt = 

(xi — i) + • • • + (x m — ~) + ~ — mt. -F satisfies the hypothesis of 
Proposition 14.31 if and only if t < 1/re. We obtain: 

(25) If cos 2 9 = t<-, ICI < *"* . 
y ' n ' ' ~ l/n -t 

4.4. Christoffel-Darboux formula and an explicit bound. It remains 
to apply the standard method with Christoffel-Darboux formula. For k = 
(ki, . . . , k m ) and I = . . . , l m ), the notation l<k stands for: l{ < ki for 
all l < i < m. 

Proposition 4.4. Let y = (yi, ... , y m ) £ R m and k = (h,... , k m ) £ N m , 
and define 

m kj 

K k (x,y) :=J2dlPl( x ) P L(v) = H (£ ^foWifc)) 

l<k j=l i=0 

and 

m kj 

N k (x,y) := ^2d kt a kt Q kt (x t ,y t )Y[{^2diPi(xj)Pi(yj)) 

t=l j^t i=0 

where 

Qk t (x t ,yt) ■= P kt +i(xt)Pk t (yt) - Pk t {x t )Pk t +i{yt)- 

Then we have the Christoffel-Darboux type formula: 

N k {x,y) 

K k {x,y) = —= — . 

o-(x) - a{y) 
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Proof. Since a{x) - a{y) = Y%L\ x t ~ Et=i Vt = YJtU( x t ~ Vt), 

m m kj 

(a(x)-a(y))K k (x,y) = - y t )) H (X> P ^M%)) 

t=l 3=1 i=0 

m kt k j 

= E - vt^diPii^Piivt)) II (E^fe)^^)) 

t=l i=0 3'^t i=0 

The Christoffel-Darboux formula for the polynomials Pk gives: 

kt 

(x t - yt)^2diPi{x t )Pi{yt) = d kt a kt Q kt (x t ,y t ) 

i=0 

with the notations of the proposition, hence the result. 

□ 

Following the standard method, we apply Proposition 14.31 to the function 

Nk(x,y) 2 
a(x) - a(y)' 

Proposition 4.5. Let y = (yi,..., y m ) G R m and k G N m ; and define 

F(x) := ^L^lE. = (o-^) _ a(y))K k (x,y) 2 = K k (x,y)N k (x,y). 

If y satisfies the conditions: 

(i) Pi{yt) > for all < i < kt and for all 1 < t < m 

(ii) P kt+ i(y t ) < for alll<t<m 

then F satisfies the hypothesis of Proposition \4-3\ for all 9 such that mt < 
cr(y). Consequently, for any code C in X m with minimum angular distance 



(26) 

(YZi d kt a kt {P kt {yt) - P kt+ i(yt)) (E?=o^(%))V 

\c\ < 



(m - a{y)) ET=i d kt a kt P kt {y t )P kt+1 {y t ) \[^ t MPiiVj)) 



Proof. Clearly, under the assumptions (i) and (ii), K k (x,y) and N k (x,y) 
have non-negative coefficients on the Pi. This is enough to ensure that it is 
also the case for the product K k (x, y)N k {x, y) (recall that the product of two 
polynomials with non-negative coefficients on the P k also has non-negative 
coefficients on the P k . This property transfers straightforwardly to the P k ; 
it is anyway general to any family of zonal polynomials). 

Obviously the sign of F(x) is the sign of cr(x) — o~(y) so the conditions of 
Proposition 14.31 are fulfilled. 
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It remains to compute /q = [F, 1] and F(l, . . . , 1). 
[F,l] = [K,iV] 

m kj m kj 

= [II ( 53 ^i-^C^J*)-^ (%*)) J 5^ dfctOfctQA^ 2/t) IX C52 d i P i( x j) P i(yj)) 
j=l i=0 t=l j^t i=0 

m m kj kj 

t=l j=l i=0 j^t i=0 

m kt 

= ^d kt a kt \^2d i P i {x t )P i (y t ),Q kt (xt,yt) 
t=l i=0 

. n [ e ^ )^ (%•)> e ^ ^ (f j ) 

j'^t i=0 i=0 

= ^4 t a fct ( - F fet (y t )P fet+ i(y t )) ]X (J^diiPiiyj)) 2 ) 

t=l j^t i=0 

where the last equality follows from (|22[). 
Let us now compute . . . , 1). We have: 

/?(1 ,..., 1) = M 

m - cr(y) 

and 

m kj 

Nkihy) = Y, d ^k t {P kt {yt) - p k t+ x(yt)) II (X> p <(f;))- 

t=i jyt i=o 

Applying the resulting bound of Proposition 14.31 leads to the announced 
bound. □ 

We proceed to choose the parameters y and k such that the conditions 
of Proposition 14.51 will be satisfied. We follow the standard method. We 
first choose the multi- index k such that rat < Y2tLi z k t i where z kt is the 
largest zero of P kt . The interlacing property of the zeros of the orthogonal 
polynomials P k guarantees that there exists y such that z kt < y t < z kt+ i 
and 

P kt (yt) + p k t +i(yt) = o 

for all 1 < t < in. Thus, Pi{yt) > and P kt +\{yt) < for all < i < k t and 
1< t < m. 
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Now we have: 

m kj 

fo = [F,l]=Y J d kt a kt {Pk t {yt)fX{{Y, d ^ P ^)) 2 ) 

t=l j^t i=0 

m 



t=l l<k 
h=k t 



and 



m kj 

Nk(L y) = 2 E d kt a kt Pk t (yt) n ( E d i p i(vj)) 

t=l j^t i=0 



rn 



2 E°^ E d L p i^- 



t=l l<k 

h=kt 

With Cauchy- Schwartz inequality (applied twice), 

E E (/ ' 

t=l Kfc 
Zt=fct 
m m 

! 

< , 

t=l Kfc 
k=k t 



< 



(m 


- *(v)) 




4 


(m 


- *(*)) 




4 


(m 


- *(*)) 




4 


(m 


-*(*)) 




4 


(m 


-*(*)) 



m m 

E^(EME rf i)(E^(y)) 2 )) 



t=l Kfc Kfc 

m 



t=l Kfc t=l Kfc 

m m kt 

E^)n(E^)/o 

t=l (=i i=o 



Denote D^ t := X^lo^- We obtain 

4(5Xia*)n£iA* 



|C|< 



m-a(y) 

We summarize the above result in the following statement: 

Proposition 4.6. For any code C in X m with minimum angular distance 6, 
for any multi-index k such that mt < J2t=i z kt> ^ Vt satisfy z kt < yt < Zfc t +i 
and 

P fet (yt) + P fct+ i(yt) = 0, 
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then 



m-a(y) 

Remark 4.7. Using the so-called adjacent polynomials instead of the Gegen- 
bauer polynomials, an enhancement of (|27|) was derived for m = 1 ■ 
It seems that this can be generalized for all m. 

4.5. A bound for the asymptotic rate. Now we consider the limit when 
n — > +oo of the rate R(C) := ^ In |C| (of course the space P 2 (0) is not 
concerned anymore) of the codes C of X m . We derive an upper bound for 
this limit from (|27|) , The next theorem settles the result obtained that way 
only in the case X = P n_1 (iC) because this bound, in the case of X = S 71 ^ 1 , 
turns out to be the same as the one obtained from the trivial isometric 
embedding (S" 1 " 1 )™ -> S mn ~ l (see Remark Oil. 

Theorem 4.8. Let C be a code in X m , X = P™" 1 ^), with minimum 
angular distance 0, and let (9i, . . . ,0 m ) £ [0, 7r/2] m satisfy Ylu=i cos2 ®t = 
mcos 2 9. Then, when n — ► oo, 

(28) R(C) < c(R LP {e x ) + ... + R LP (e m )), 

where Rlp is defined in (|2.2[l . 

Proof. Same as in J3], involving the asymptotic estimate of Zf.. We reproduce 
it here: Consider an infinite sequence k(n) such that 2k(n)/cn tends to a 
finite limit p as n tends to infinity. Then 01 

t bmz fc(n) =4 ^_ 1 + 2)2 



and, since from jS], 



we have 



0l = ,'f»+* 



1 2 f-n + kin) — 1 \ 

lim -hxD k t n) = lim -In 2 // = c((l + p) ln(l + p) - plnp). 

n^oo n n— >oo n y kin) I 



Inverting the conditions 



(Pt~ 1 + 2) 2 
leads to 

1 — sin t 
2 sin 0t 

Let fe( = [p t n\ , and let y t satisfy < y t < z kt+ i and i\ (y t )+Pk t +i(yt) = 
(the existence of yt is guaranteed by the interlacing property of the zeros 
of the Jacobi polynomials). Then from (|27|). 

, c , < 4 (IXi ak t )UT=i D k t 
m-a(y) 
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Since o~(y) ~ m cos 2 9 and the expression has a finite limit when 

kf/n tends to pt, the rate -R(C) satisfies 

1 m m 



n 

t=i t=l 



□ 



Remark 4.9. • It is worth noticing that the choice 9t = in (j28|) 

yields to the bound 

(29) i?(C) < cmR L p(9). 

This bound can be derived more easily, since every code in X m is 
also a code in the cmn-th dimensional unit sphere (combining the 
mapping f3 form = 1 and the obvious mapping (S" 1-1 ) — > S mn ^ 1 ). 
It turns out that, since the function Rlp(0) as a function of t = 
cos 2 9 is not convex, the bound ({2"%)) slightly improves on We 
discuss this in more details in the next subsection. 
• The same method applied to X = S^ 1 would lead to: 

m 

R(C) < Rlp{9\) + . . . + RLp{9m), f or oil 9 t such that ^^cos6>j = mcos9. 

t=i 

But the function Rlp{9) as a function of t = cos9 is convex, there- 
fore the choice of (9\, . . . , 9 m ) that minimizes the right hand side is 
6»i = • • • = 9 m = 9, yielding (gHJ). 

4.6. Analysis of ((2"5|) versus (|29j) . Let C 2 be the set of continuous, twice 
differ entiable functions with continuous second derivative. For a function / 
defined on [0, 1[, of class C 2 , we denote: 



f (m \t) := min 



ti,...,t m e[o,i[ m 

Clearly, if / is convex on [0, 1[, we have f( m ^ = f, and, if / < g, f^ < g( m \ 
It is also easy to see that f( m ' < /( m ) when m divides m'. 

The function we are interested in is f(t) = Rlp{9) where t = cos 2 9. We 
have 

f(t) = (1 + p{t)) ln(l + p{t)) - p{t) ln(p(t)) 

where 

P{t) = \(-l + (l-t)-^ 2 ). 

One can check that the second derivative of / takes negative values on 
some interval [0, to], to — 0.208, and then takes positive values on [to>l[- 
The function / is an increasing function, with /(0) = 0, first concave then 
convex. We consider the function g on [0, 1[, whose graph C g determines the 
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convex hull of the portion of plane above the graph Cf of /. The function g 
is uniquely determined by the conditions: 

9<f 
g is convex 

g is maximal with these properties 

Let us denote by t\ the unique value for which the tangent at (t%, f(t\)) to 
Cf contains the origin (0, 0). The value t\ ~ 0.379 is the unique solution to 

/(*) = f(t)t 

and the slope of the tangent to Cf at t\ equals f'(ti) ~ 1.089. Then the 
function g is defined by: 

g (t) = f'it^t ~ 1.089* for all t £ [0,h] 
g(t)=f(t) for all t G [t±, 1[ 

Since g is convex, we have for all m and all t £ [0, 1[, g^ m \t) = g(t) < f^ m '{t). 
In other words, on [0, ti], /( m ) is somewhere between g and /, and on [t±, 1[, 
j(m) _ j _ Clearly, when m — > +oo, 

f(m) _^ g Als0; the 

maximum 5 of 

f(t) — g(t) is an upper bound for the maximum of f(t) — y( m ) (f ). Numerical 
calculation gives 5 ~ 0.016. Considering our primary goal, i.e., to compare 
((25|) and this means that the improvement of upon is upper- 
bounded by 0.016m. 

It seems difficult to determine the optimal choice of (ti,...,t m ) that 
minimizes the quotient ^ tl ^ + m + ^ fm ^ . A natural choice is (ti, . . . ,t m ) = 
(0, 0, ..,mt/r, . . . , mt/r) with r non-zero and equal coordinates. In that case, 

m)+ m +/( * m) = kf^f) and re q uires * < r l m - If * = it is certainly the 
best choice since then the resulting point lies on C g . Numerical experiments 
seem to show that, for m = 2,3, and t < l/m, r = 1 does minimize the 
quotient /Mt^) , 

5. Bounds for codes in the Grassmann and Stiefel manifolds 

In this section, we summarize the consequences of the above results for 
Grassmann and Stiefel codes. Following a standard notation in coding the- 
ory, we denote by A(X, d), the maximal number of elements of a code C of 
the space X with minimum distance d. 

We have proved in the section 13.31 that the size of Grassmannian codes 
with minimal chordal distance d = \Jm — s is upper bounded by the size of 
codes in ¥ n ~ l (K) m with minimal angular distance 9, where cos 2 6 = s/m. 
Thus we have proved that: 



(30) A(g m , n (K), d) < A(P n - 1 (K) m , 9) with 9 = arccos y / l-d 2 /m. 

Linear programming bounds on AiW"""" 1 (K) m , 9) were derived in Section 0J 
We believe that these bounds are not good in general for finite val- 
ues of the parameters, because we use only a rough estimate of a(p, q) = 
trace(7T p o ir q ) in the inequality (|TH)) (we replace Si<ij<m( e * ' e 'j) 2 wr th 
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J2i<i< m ( e i ' e 'i) 2 )- -^ we compare the bounds obtained with the zonal poly- 
nomials of small degree, (|25|) is worse than the simplex bound, obtained 
from the zonal polynomial of degree 1 of G m ,n(M)- Moreover, numerical ex- 
periments for small parameters m and n (with the package LRS, by David 
Avis, http://cgm.cs.mcgill.ca/~avis/C/lrs.html), confirms that the bounds 
obtained from the zonal polynomials of G m! n(M) are sharper than the ones 
obtained from Proposition Ol for X = P™- 1 ^)" 1 . 

Surprisingly, the consideration of P n_1 (i^) m allows us to obtain better 
bounds for the asymptotic rate than the ones obtained previously by either 
the isometric embedding given in |21| of G m ,n into a unit sphere of the 
dimension (n — 1) (n + 2)/2 (see also Remark 13.3(1 . or the spectral method 
developed in |Hj with the zonal polynomials of Q m n . We summarize the new 
bound we have obtained in the next theorem: 

Theorem 5.1. Let C be a code iti Gm,n(,^-') with minimal chovdol distance 
d = y/m — s, and let 9 = arccos(y / s/m). Then, when n — > +00, 

(31) R(C) <mm{R 1 (d),R 2 (d)}, 
where 

(32) Ri(d)= min c(R LP (e 1 ) + --- + R L p(e m )) 
and 

(33) R 2 (d) = cmR s {9). 

The bounds R\{d) and R 2 (d) are depicted in Figure [21 
We have proved in Section 13.31 that the size of Stiefel codes with minimal 
chordal distance d = \f2\Jm — s is upper-bounded by the size of codes in 
(S"^ -1 ) with minimal angular distance 6, where cos 9 = s/m. In Section0J 
we derived linear programming bounds for codes in ('5 ,cn - 1 ) m ) thus implying 
bounds for Stiefel codes. These bounds are, up to our knowledge, the first 
general bounds for Stiefel codes, and we believe that they are rather sharp. 
For the asymptotic rate, the best result is obtained in (|2Uj) . We summarize 
these results in the following theorem: 

Theorem 5.2. With the previous notations, and 9 = arccos (1 - d 2 /2m), 

(i) A(V m>n (K),d) < A{{S™~ l ) m ,9) 

(ii) When n -» +oo ; R{C) < cmR s {9) 

6. Conclusions 

Using relations between Grassmann and Stiefel manifolds and other spa- 
ces, we derive new bounds on the size of Grassmannian codes (Theorem 
15. 1|) and Stiefel codes (Theorem 15 .2(1 . These are the best known asymptotic 
bounds on the rate of Grassmannian and Stiefel codes. 
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1.1 1.2 1.3 1.4 1.5 1.6 1.7 



Figure 2. Upper bounds on the asymptotic rate of real 
Grassmannian codes with minimum chordal distance d, m = 
3. The solid line is Ri(d), and the dashed line is R,2(d) (see 
Theorem E2J) 
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